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Abstract 

We derive a simple relation for estimating the relative emittance growth in x and y due to 
intrabeam scattering (IBS) in electron storage rings. We show that IBS calculations for the ATF 
damping ring, when using the formalism of Bjorken-Mtingwa, a modified formalism of Piwinski 
(where n 2 / (3 has been replaced by TC), or a simple high-energy approximate formula all give results 
that agree well. Comparing theory, including the effect of potential well bunch lengthening, with 
a complete set of ATF steady-state beam size vs. current measurements we find reasonably good 
agreement for energy spread and horizontal emittance. The measured vertical emittance, however, 
is larger than theory in both offset (zero current emittance) and slope (emittance change with 
current). The slope error indicates measurement error and/or additional current-dependent physics 
at the ATF; the offset error, that the assumed Coulomb log is correct to within a factor of 1.75. 
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I. INTRODUCTION 



In future e + e~ linear colliders, such as the JLC/NLC0, §, damping rings are needed 
to generate beams of intense bunches with low emittances. The Accelerator Test Facility 
(ATF) || at KEK is a prototype of such damping rings. One of its main goals, and one that 
has been achieved, was the demonstration of extremely low vertical emittances^, 0]. At the 
low ATF emittances, however, it is found that intrabeam scattering (IBS) is a strong effect, 
and one that needs to be understood. 

Intrabeam scattering is an effect that depends on the ring lattice — including the errors — 
and on all dimensions of the beam, including the energy spread. At the ATF all these 
dimensions can be measured; unique to the ATF is that the beam energy spread, an espe- 
cially important parameter in IBS theory, can be measured to an accuracy of a few percent. 
In April 2000 the single bunch energy spread, bunch length, and horizontal and vertical emit- 
tances were all measured as functions of current over a short period of time|6|, [?J. The short 
period of time was important to ensure that the machine conditions remained unchanged; the 
bunch length measurement was important since potential well bunch lengthening is signifi- 
cant at the ATF |7j . The question that we attempt to answer here is, Are these measurement 
results in accord with IBS theory? 

Intrabeam scattering theory was first developed for accelerators by Piwinski[^[, a result 
that was extended by Martini||, to give a formulation that we call here the standard Piwin- 
ski (P) method[T(J; this was followed by the equally detailed Bjorken and Mtingwa (B-M) 
result ||TT|. Both approaches solve the local, two-particle Coulomb scattering problem for 
(six-dimensional) Gaussian, uncoupled beams, but the two results appear to be different; 



of the two, the B-M result is thought to be the more general [0. Other simpler, more 
approximate formulations developed over the years are ones due to Parzen|l3|, Le Duff|jl4 



Raubenheimer |15| , and Wei]T6[. Recent reports on IBS theory include one by Kubo and 
Oide, who adapt an intermediate result from Bjorken- Mtingwa's paper to find the solution 



for cases of arbitrary coupling [17]], a method that is now used in the optics computer pro- 
gram SADjTBj; and one by Venturini that solves for IBS in the presence of a strong ring 
impedance |1^ . 

Intrabeam scattering measurements have been performed primarily on hadronic ][20|, |21| , 
22 , p3| and heavy ion machines p| , p5|, where the effect tends to be more pronounced, 
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though measurement reports on low emittance electron rings can also be found |26|, |27 . 
Typical of such reports, however, is that although good agreement may be found in some 
beam dimension(s), the set of measurements and/or agreement is not complete (e.g. in 
Ref. [M growth rates agree reasonably well in the longitudinal and horizontal directions, 
but completely disagree in the vertical). Note that one advantage of studying IBS using 
electron machines is that it can be done by measuring steady-state beam sizes. 

In this report we briefly describe intrabeam scattering formulations, apply and compare 
them for ATF parameters, and finally compare calculations with the full set of data of April 
2000. For more details on the hardware and such measurements at the ATF, the reader is 
referred to Ref. @, |[]. 

II. IBS CALCULATIONS 



We begin by describing the method of calculating the effect of IBS in a storage ring. Let 
us first assume that there is no x-y coupling. 

Let us consider the IBS growth rates in energy p, in the horizontal x, and in the vertical 
y to be defined as 

J_ l_dop 1_ J_(k^_ 1_ 1 dt l J 2 

T p ~V p Ht' Y X ~^^T' ¥ y ~^^F- u 

Here a p is the rms (relative) energy spread, e x the horizontal emittance, and e y the vertical 
emittance. In general, the growth rates are given in both P and B-M theories in the form: 

^ = </<> (2) 

where subscript % stands for p, x, or y. The functions fi are integrals that depend on 
beam parameters, such as energy and phase space density, and lattice properties, including 
dispersion; the brackets () mean that the quantity is averaged over the ring. In this report 
we will primarily use the fi of the B-M formulation ||35||. 

From the 1/Tj we obtain the steady-state properties for machines with radiation damping: 
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where subscript represents the beam property due to synchrotron radiation alone, i.e. in 
the absence of IBS, and the r, are synchrotron radiation damping times. These are 3 coupled 
equations since all 3 IBS rise times depend on e x , e y , and a p . 



The way of solving Eqs. ^ that we employ is to convert them into 3 coupled differential 
equations, such as is done in e.g. Ref. PB| , and solve for the asymptotic values. For example, 
the equation for e y becomes 

de y _ 2(e y - e y0 ) 2e y 

it- — + t; • (4) 



and there are corresponding equations for e x and a 



Before solving these equations one needs to know the source of the vertical emittance 
at zero current. We consider 3 possible sources: (i) vertical dispersion due to vertical orbit 
errors, (ii) (weak) x-y coupling due to such things as rolled quads, etc, and (iii) a combination 
of the two. If the vertical emittance at zero current is due mainly to vertical dispersion, 
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e y0 ^Je(n y )a^ , (5) 

with J t the energy damping partition number and 7i = [rf + (firf — ^f3'r)) 2 ]/f3 the dispersion 
invariant, with rj and (3, respectively, the lattice dispersion and beta functions. If e y o is 
mainly due to coupling we drop the e y differential equation and simply let e y = Ke x , with 
k the coupling factor. In case (iii) we approximate the solution by replacing the parameter 
e y o in Eq. [| by the quantity [/te x (l — r y /T y ) + e y0( i\, where e y od is the part of e y o due to 
dispersion only. Note that the practice — sometimes found in the literature — of solving IBS 
equations assuming no vertical errors, which tends to result in near or even negative 
vertical emittance growth, may describe a state that is unrealistic and unachievable. Note 
also that in case (i) once the vertical orbit — and therefore (H y ) — is set, e y o is no longer a 
free parameter. 

In addition, note that: 

• A fourth equation in our system, the relation between bunch length a s and a p , is also 
implied; generally this is taken to be the nominal (zero current) relation. In the ATF 
strong potential well bunch lengthening, though no microwave instability, is found at 
the highest single bunch currents]?]. In our comparisons with ATF measurements we 
approximate this effect by adding a multiplicative factor f pw (I) [I is current], obtained 
from measurements, to the equation relating a s to a p . (Note that potential well bunch 
lengthening also changes the longitudinal bunch shape, a less important effect that we 
will ignore.) 
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The B-M results include a so-called Coulomb log factor, \n(b max /b m i n ), with b max , 
b m i n maximum, minimum impact parameters, quantities which are not well defined. 



For round beams it seems that b max should be taken as the beam size ^9|. For bi- 
Gaussian beams it is not clear what the right choice is. Normally b max is taken to be 
the vertical beams size, though sometimes the horizontal beam size is chosen pt)|]. We 



take b max = a y ; b m i n = r oc 2 /(vl) = r^ x j{pi 1 e x \ with r the classical electron radius 
(= 2.82 x 10~ 15 m), v x the transverse velocity in the rest frame, and 7 the Lorentz 
energy factor. For the ATF, the Coulomb log, (log) = 16.0. 

The IBS bunch distributions are not Gaussian, and tail particles can be overemphasized 
in these solutions. We are interested in core sizes, which we estimate by eliminating 



interactions with collision rates less than the synchrotron radiation damping ratepl 
We can approximate this in the Coulomb log term by letting 7r&^ w - n (| , y a .|)(n) equal the 
synchrotron damping rate in the rest frame, with n the particle density in the rest 



frame fPT|; or b min = J AT[a x a y a z ^j /[Nct](/3 x / 'e a .) 1//4 , with N the bunch population. For 
the ATF with this cut, (log) = 10.0. 



A. High Energy Approximation 



For both the P and the B-M methods solving for the IBS growth rates is time consuming, 
involving, at each iteration step, a numerical integration at every lattice element. A quicker- 
to-calculate, high energy approximation, one valid in normal storage ring lattices, can be 



derived from the B-M formalism 32 
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The requirement on high energy is that a,b -C 1; if it is satisfied then the beam momentum 
in the longitudinal plane is much less than in the transverse planes. For flat beams a/6 is 



less than 1. In the ATF, for example, when e y /e x ~ 0.01, a ~ 0.01, b ~ 0.1, and a/b ~ 0.1. 
The function g, related to the elliptic integral, can be well approximated by 

g(a) « a (o.o2i-o.o44m Q ) [for Q1 < a < . (g) 

to obtain g for a > 1, note that (7(a) = g(l/a). 

Note that Parzen's high energy formula is a similar, though more approximate, result 
to that given here [13]; and Raubenheimer's approximation is formulas similar, though less 



accurate, than the first and identical to the 2nd and 3rd of Eqs. P|l5|. Note that Eqs. || 
assume that e y0 is due mainly to vertical dispersion; if it is due mainly to x-y coupling we let 
7iy = 0, drop the l/T y equation, and simply let e y = ne x . Finally, note that these equations 
still need to be iterated, as described before, to find the steady-state solutions. 

B. Emittance Growth Theorem 



Following an argument in Ref. [IS| we can obtain a relation between the expected vertical 
and horizontal emittance growth due to IBS in the presence of random vertical dispersion: 
We begin by noting that the beam momentum in the longitudinal plane is much less than 
in the transverse planes. Therefore, IBS will first heat the longitudinal plane; this, in turn, 
increases the transverse emittances through dispersion (through TC, as can be seen in the 
2nd and 3rd of Eqs. |]), like synchrotron radiation (SR) does. One difference between IBS 
and SR is that IBS increases the emittance everywhere, and SR only in bends. We can write 

tyO ^ 3x (7~ty)b £y ~ e y0 ^ Jx\%y) (10) 

£x0 <JyV~ix)b — CxO <Jy(Hx) 

where J x%y are damping partition numbers, and ()& means averaging is only done over the 

bends. For vertical dispersion due to errors we expect (H y )b ~ (Hy)- Therefore, 

r _ ( e y ~ e yo)/ e yQ _ (?Qb ^ 
(e x — £ x o)/txo (T~tx) 

which, for the ATF is 1.6. If, however, there is only x-y coupling, r t = 1; if there is both 
vertical dispersion and coupling, r e will be between (TC x )b/('H x ) and 1. 

C. Numerical Comparison 



Let us compare the results of the methods P, B-M, and Eq. |6] when applied to the ATF 
beam parameters and lattice, with vertical dispersion and no x-y coupling. We take as 
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parameters those given in Table I, and, for this comparison, let f pw = 1. In addition we 
have J e = 1.4, (f3 x ) = 3.9 m, (/3 y ) = 4.5 m, (i] x ) = 5.2 cm and (TC X ) = 2.9 mm. To generate 
vertical dispersion we randomly offset magnets by 15 /xm, and then calculate the closed orbit 
using SAD. For our seed we find that the rms dispersion (f] y ) r ms — 7.4 mm, (H y ) = 17 /im, 
and e y o = 6.9 pm (in agreement with Eq. [5]). For consistency between the methods we here 
take the cut-off parameter d = 3a y in P to corresponds to (log) = In [da^ / (Ar^a 2 )} = 16 in 



B-M. 

TABLE I: Typical ATF parameters in single bunch mode. 



Circumference 


C 


138 


m 


Energy 


E 


1.28 


GeV 


Current 


I 


3.1 


mA 


Nominal energy spread 


CpO 


5.44 


io- 4 


Nominal horizontal emittance 


CxO 


1.05 


nm 


Nominal bunch length 


&s0 


5.06" 


mm 


Longitudinal damping time 


Tp 


20.9 


ms 


Horizontal damping time 


T x 


18.2 


ms 


Vertical damping time 


T V 


29.2 


ms 



a at rf voltage 300 kV 



Performing the calculations, but first comparing the standard Piwinski and the B-M 
methods, we find that the growth rates in p and x agree well; the vertical rate, however, 
does not. In Fig. |I| we display the local IBS growth rate in y over half the ring (the periodicity 
is 2), as obtained by the two methods, and see that the P result, on average, is 25% low. 
Studying the two methods we note that a conspicuous difference between them is their 
dependence on dispersion: for P the fi depend on it only through r] 2 /f3; for B-M, through 
(ft = [rj' — ^P'rj/fl] and through H. Let us replace if / (5 in P with TC to create a method that 
we call the modified Piwinski result. In Ref. |$2| it is shown that, in a normal storage ring 



lattice, at high energies, the results of this method become equal to those of B-M. 

Comparing with this method we find that, indeed, the three growth rates now agree 
reasonably well with the B-M result. Fig. [| displays the 3 local growth rates as obtained by 
the modified P and B-M methods. The 1/Tj, the average values of these functions, are given 
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FIG. 1: Vertical steady-state (local) growth rate over 1/2 the ATF for an example with vertical 
dispersion due to random errors. Given are results due to standard Piwinski (green) and Bjorken- 
Mtingwa (blue). 

in Table II. We note that the P results are all slightly low, by 4.5%. The B-M method gives: 
Cp/cpo = 1-52, e x /e x o = 1.90, e y /e y o = 2.30. Note that for this error seed the emittance 
growth ratio of Eq. |TT] is r e = 1.44, close to the 1.6 expected for the ATF lattice. 
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FIG. 2: Steady-state (local) growth rates over 1/2 the ATF for an example with vertical dispersion 
due to random errors. Given are results due to modified Piwinski, Bjorken-Mtingwa, and Eqs. @. 



Repeating the calculation using Eqs. |6| we find that the computing time is greatly reduced, 
and the growth rates agree quite well with the B-M results (see Table II). The dots in Fig. |2| 
give the local rates corresponding to Eqs. ||, and we see that even these agree quite well. 
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TABLE II: Steady-state IBS growth rates for an example including vertical dispersion due to 
random errors. 



Method 1/T P [s- 1 ] 1/T X [s- 1 ] 1/T y [s" 1 ] 

Modified Piwinski 25.9 24.7 18.5 

Bjorken-Mtingwa 27.0 26.0 19.4 

Eqs. | 27.4 26.0 19.4 



D. Comparison with SAD Results 

The optics program SAD basically follows the B-M formalism, but it does it in a form 
that treats the three beam directions on equal footing. The final results are given in terms of 
the normal modes of the system and not the beta and dispersion functions of the uncoupled 
system (as in our approximation). For vertical dispersion dominated problems there is no 
difference in result. In coupling dominated problems there will be a difference that, in the 
case of small x-y coupling due to errors, we expect to be small. 

We consider the ATF lattice with random magnet offsets and rotations. Other machine 
parameters are the same as before; again / = 3.1 mA. For this lattice (f] y ) r m,s = 7.4 mm 
and e y0 /e x0 = 1%. For this problem we solve IBS using SAD (for 2 different seeds), and 
also our approximate method where we include vertical dispersion (as before) and a global 
coupling parameter k = 0.34%. We take (log) = 9.1. Comparing steady-state local growth 
rates, we find good agreement in p and x for all three calculations. In y, however, there is 
a significant variation (see Fig. ^j. The growth rates, the average values of these functions, 
however, agree well (see Table III). Note that the steady-state relative growths in (a p ,e x ,e y ) 
are (1.38,1.56,1.64) for SAD, and (1.38,1.62,1.61) for our approximate calculation. 

TABLE III: Steady-state IBS growth rates for an example including vertical dispersion and x-y 
coupling due to random errors. 

Method l/T p [s- 1 ] 1/T X [s- 1 ] 1/T y [s" 1 ] 

SAD, seed 1 22.5 19.6 13.1 

SAD, seed 2 22.3 19.6 13.5 

Our approx. calculation 22.9 21.0 12.9 
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60 

SAD, seed 1 

40 

SAD, seed 2 
■ - ■ coupling approx.^ 



20 40 60 

s [m] 

FIG. 3: Vertical steady-state (local) growth rate over 1/2 the ATF for an example with vertical 
dispersion and x-y coupling due to random errors. Given are results obtained by SAD (for 2 seeds; 
the solid curves) and by the coupling approximation used here (the dots). 

III. COMPARISON WITH MEASUREMENT 
A. Measurements 

At the ATF the energy spread and all beam sizes can be measured. Unique at the ATF is 
that the energy spread, a particularly important parameter in IBS theory, can be obtained 
to a few percent accuracy. In this measurement the beam is extracted and its size measured 
on a screen in a highly dispersive region. The bunch length is determined with a streak 
camera in the ring. The emittances can be measured using 3 methods: wire monitors in the 
extraction line, a laser wire in the ring, and an interferometer in the ring. Unfortunately, 
for e y all 3 methods have their difficulties. For example, the wire measurement is very 
sensitive to optics errors (such as roll and dispersion) in the extraction line. Or, the laser 
wire measurement, being time consuming (taking ~ 1 hour per measurement), is sensitive 
to drifts in machine and beam properties. 

Because of the effects of IBS the energy spread measurement (which is quick and easy 
to perform) has become a useful technique for monitoring changes in beam size. Thus, 
evidence that we are truly seeing IBS at the ATF include: (1) when moving onto the coupling 
resonance, the normally large energy spread growth with current becomes negligibly small; 
(2) if we decrease the vertical emittance using dispersion correction, the energy spread 
increases. 
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B. Comparison with Theory 



In Fig. as an example, we present the time development, after injection, of energy 
spread for 3 different beam currents (the plotting symbols). The measurement was per- 
formed by continually injecting beam into the ATF, while varying the extraction timing. 
If we take the B-M formalism, with f pw = 1, and with x-y coupling 0.006, and solve the 
differential equations for energy spread and beam sizes, we obtain the curves in the fig- 
ure (if we include potential well distortion the fitted coupling becomes 0.0045). The short 
time (^S 0.05 s) behavior does not agree with the data, since the beam in reality enters 
the ring badly mismatched (a region which would be difficult to simulate); in the longer 
time range, however, after > 3r p , the agreement becomes quite good. The minimum in the 
curves can be explained as follows: Initially the energy spread and beam sizes reduce due to 
synchrotron radiation; when the beam volume becomes smaller than a certain amount, the 
energy spread begins to increase due to IBS. This result indicates reasonably good agreement 
between measurement and theory. 
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FIG. 4: Measured energy spread as function of time after injection, for 3 different currents (the 
plotting symbols). The curves give B-M simulations assuming an x-y coupling of 0.006 and no 
potential well distortion. This plot is reproduced from Ref. 0]. 

To compare with theory absolutely, however, we need to measure all beam properties with 
the machine in the same condition. Such a complete series of measurements was performed 
on stored beam at the ATF over a short period of time in April 2000. The rf voltage was 
V c = 300 kV. The energy spread and bunch length vs. current measurements are shown in 
Fig. m The curves in the plots are fits that give the expected zero current result. Emittances 
were measured on the wire monitors in the extraction line (the symbols in Fig. P>c; note 
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that the symbols in Fig. |]a reproduce the fits to the data of Fig. ||). We see large growth also 
in the emittances. Unfortunately, we have no error bars for the emittance measurements, 
though we expect the random component of errors in y to be 5-10%, and less in x. Note 
that e y0 appears to be about 1.0-1.2% of e x0 . 

B 

ST 6 
o 

^ 4 
b 

2 


0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 
I [mA] I [mA] 

FIG. 5: Measurements of steady-state energy spread (a) and bunch length (b), with V c = 300 kV. 
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0._ _ _ ... 

I [mA] 



I [mA] 

FIG. 6: ATF measurement data (symbols) and IBS theory fits (the curves). The symbols in (a) 
give the smooth curve fits to the measured data of Fig. |5[ 

Let us compare B-M calculations with the data. Here we take f pw as given by the 
measurements, and take as Coulomb log our best estimate, (log) = 10. Note that in the 
machine the residual dispersion is typically (r] y ) rms ~ 3 mm. To set our one free parameter, 
e y0 , we adjust it until at high current a p agrees with the measurement. In Fig. ^] we give 
examples: 

1. Vertical dispersion only, with (r} y ) rms = 5.6 mm and e y o = 4.0 pm (solid); 

2. Coupling dominated with k = 0.33% (dashes); 
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Example (1) 
Example (2) 
Example (3) (IBS 
stronger by 1.75) 
Ex. (2) + 0.9% x-y 




coupling meas. error n.O 0.5 1.0 1.5 2.0 2.5 3.0 



3. Coupling dominated with k = 1.2%, with the Coulomb log artificially increased by a 
factor 1.75 (dotdash); 

4. Same as Ex. 2 but assuming e y measurement error, i.e. adding 0.9% of the measured 
(and splined) e x to the calculated e y (the dots). 

We see that cr p (I) agrees well with the measurements for all cases, and e x (I) agrees reason- 
ably well. For Examples 1 and 2, however, e y o is significantly lower than the measurements 
seem to indicate, and the growth with current is also less. To obtain reasonable agreement 
for e y0 we need to assume that either IBS is ~ 75% stronger (in growth rates) than theory 
predicts, or there is significant measurement error, equivalent to ~ 1% x emittance coupling 
into the y measurement. Yet even with such assumptions the e y (I) dependence does not 
agree. 

What does the emittance growth theorem of Sec. |IIB| say about these results? It appears 



that e x grows by ~ 85% by / = 3 mA; e y begins at about 1.0-1.2% of e x o, and then grows to 
about 3% of e X Q. Therefore, the relative emittance growth ratio is r t ~ 2.1-2 A, much larger 
than the expected result if we are coupling dominated (1.0); and still significantly larger 
than the expected result if we are dispersion dominated (1.6), a case that is anyway unlikely 
since it requires an implausibly large (r] y ) rms ~ 9 mm. Thus, the emittance growth theorem 
indicates that e y (I) as measured is not in agreement with IBS theory. 



IV. DISCUSSION 



Our disagreement in e y between theory and measurement consists of two parts, an offset 
part (e y o) and a disagreement in slope (de y /dl). Together they indicate that we have: error 
in theory, additional physics at the ATF, and /or error in measurement. 

IBS theory is a mature theory, and the relation between longitudinal and transverse 
growth rates (the 2nd and 3rd of Eqs. []) is simple and intuitively easy to understand. The 
main uncertainty in theory may be with the scale factor, particularly in the Coulomb log 
factor for beams with elliptical cross-section. Yet a scale factor error can affect only the 
offset part of the disagreement. Note also that even if the argument of (log) were in error 
by an order of magnitude this part of the disagreement would be changed by only a small 
amount (25%). 
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The disagreement in de y /dl might be explained by the presence of additional current- 
dependent physics at the ATF. We have seen that <r p (I) and e x (I) can be made to agree 
reasonably well between theory and measurement; at the same time, however, the mea- 
sured e y (I) grows much faster than predicted. One might, therefore, suspect the presence 
at the ATF of another current-dependent effect, one that increases the projected vertical 
emittance — though not the real emittance. An example of such an effect is a y-z tilt of the 
beam induced by closed orbit distortion in the presence of a transverse impedance [|l|, |33 |. 
More study needs to be done in this direction. 

As mentioned before, measuring accurately the small vertical emittances at the ATF is 
difficult, and, therefore, emittance measurement error is likely responsible for much of the 
disagreement found. We noted that a 1% coupling measurement error in the extraction line 
wire measurements can account for the offset part of the disagreement; the slope disagree- 
ment, however, is not easy to explain assuming measurement error alone (for an attempt in 
this direction, see e.g. Ref. [Q). 

Over the time since April 2000 the systematics of the emittance measurements have 
improved, especially for the laser wire measurement. Newer results seem to suggest that the 
April 2000 measured vertical emittance may have been too large 0, |5|. For the near future 
we urge that the effort to obtain reliable emittance measurements at the ATF be continued. 
In addition, experiments to study the possible existence of other current-dependent effects 
should also be performed. Ultimately, one goal should be to test the accuracy of theoretical 
IBS growth rates to the 10-20% level. Note that once we are successful at such benchmarking 
experiments, we will be able to use the ATF energy spread measurement as a diagnostic for 
the absolute emittances of the beam. 



V. CONCLUSION 



We began by describing intrabeam scattering calculations for electron storage rings, fo- 
cusing on machines with small random magnet offset and roll errors. We derived a simple 
relation for estimating the relative emittance growth in x and y due to IBS in such machines. 
We have shown that IBS calculations for the ATF damping ring, when using the formalism 
of Bjorken-Mtingwa, a modified formalism of Piwinski (where rf / {3 has been replaced by Ti), 
or a simple high-energy approximate formula all give results that agree well. By comparing 
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with numerical results from SAD we have demonstrated that weak coupling due to random 
magnet roll can be approximated by solving the uncoupled problem with the addition of a 
global coupling parameter. 

Comparing the B-M calculations, and including the effect of potential well bunch length- 
ening, with a complete set of ATF steady-state energy spread and beam size vs. current 
measurements we have found reasonably good agreement in energy spread and horizontal 
emittance. At the same time, however, we find that the measured vertical emittance is 
larger than theory in both offset (zero current emittance) and slope (emittance change with 
current). The slope error indicates measurement error and/or the presence of additional 
current-dependent physics at the ATF. The offset error suggests that the assumed Coulomb 
log is correct to within a factor of 1.75 (though we believe that it is, in fact, more accurate, 
with part of the discrepancy due to measurement error). More study is needed. 
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